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Abstract 

The complete spectrum of collective modes of the triplet order parameter in the superfluid 
neutron matter is examined in the BCS approximation below the pair-breaking threshold. The 
dispersion equations both for the unitary and nonunitary excitations are derived and solved in the 
limit of q — > by taking into account the anisotropy of the energy gap for the case of P-wave 
pairing. By our analysis, there is only one Goldstone mode which is associated with the broken 
gauge symmetry. We found no additional Goldstone modes associated with the broken rotational 
symmetry but found that the oscillations of the total angular momentum are qualitatively similar 
to the "normal-flapping" mode in the A- phase of superfluid Helium. There are also two collective 
modes associated with internal vibrations of the structure of the order parameter oscillating with 
u (T = 0) = 1.20A and u (T = 0) = 0.61A . 
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I. INTRODUCTION 



The neutron star evolution is governed mostly by microscopic processes occuring in its 
volume and depend substantially on the spectrum of thermal excitations that can exist in 
the bulk baryon matter. Strong and electromagnetic interactions form numerous collective 
excited states in the normal (nonsuperfluid) component of the neutron star matter which 
have been extensively studied by many authors. 12]. 

The purpose of this paper is to present the complete spectrum of collective modes in 
the superfluid phase of superdense neutron matter. It is generally accepted that the pair 
condensation in superdense neutron matter occurs into the 3 P 2 state (with a small admixture 
of 3 i<2 ). Such a model is based on the properties of the bare NN interaction, which contains 
a relatively strong spin-orbit interaction in this channel. Strictly speaking, it is not really 
clear that a 3 P 2 pairing instability will survive at the relevant densities in neutron stars 
(double nuclear saturation density, or more), due to uncertainties in the 3N interaction. 
Nevertheless this model is conventionally used for estimates of neutrino energy losses in the 
minimal cooling scenarios of neutron stars. 

We do not consider the well discussed processes of pair breaking and formation [l3j, 
but focus on the undamped collective oscillations of the order parameter. The collective 
excitations below the pair-breaking threshold can play an important role in the thermo- 
dynamic properties of the neutron liquid at lowest temperatures when the other possible 
excitations are strongly suppressed by the superfluidity. Some properties of the anisotropic 



Goldstone mode in the neutron condensate were discussed in Ref. jl5[], and the existence 
of two spin- wave modes was predicted recently [l^-[l9| in the average-angle approximation. 
In the present paper we consider the solution to the dispersion equations by taking into 
account the anisotropy of the energy gap and find the whole spectrum of the undamped 
collective oscillations of the 3 P 2 order parameter. 

I — I 

In applications to neutron stars it is customary [20[ to consider the 3 Po state with a 



preferred magnetic quantum number M = 0. Sophisticated calculations [2l|-[23| have shown 
that, besides the above one-component state, there are also multicomponent states involving 
several magnetic quantum numbers that compete in energy and represent various phase 
states of the condensate dependent on the temperature. A delicate difference in the above 
gap magnitudes occurs owing to small tensor interactions between paired quasiparticles. If 



we neglect the tensor interactions the multicomponent states are degenerate and can be 
obtained from the state M = by trivial rotation of the frame in space. We therefore will 
focus on the 3 P2 neutron pairing into a state with M = 0. 

Previously eigen-modes of the order parameter have been thoroughly studied in the su- 



perfluid liquid 3 He 



24]- 



30]. The complete analysis of the order-parameter collective modes 
in anisotropic 3 He-A is given in Ref. 3jJ. The pairing interaction in 3 He is invariant with 
respect to the rotation of spin and orbital coordinates separately. This admits spin fluctua- 
tions independent of the orbital coordinates. In contrast, the spin-triplet neutron condensate 
arises in the high-density neutron matter mostly owing to the attractive spin-orbit interac- 
tions which do not possess the above symmetry. Therefore the results obtained for liquid 
3 He can not be applied directly to the spin-triplet neutron superfluid, where the most at- 
tractive channel of interactions corresponds to spin, orbital, and total angular momenta 
S = 1, L = 1, and J = 2, respectively, and pairs quasiparticles into the 3 P2 states with the 
projection of the total angular momentum M — 0, ±1, ±2. 

As is well known the strong interactions in the neutron matter are not restricted by 
the pairing forces. One can expect that the Fermi-liquid effects are able to renormalize 
the collective modes. We hope however that the principal effect occurs in the case of the 
Goldstone mode, where the sound velocity is strongly renormalized by the particle-hole 



151 ] . Since the Landau parameters 



interactions. This effect was already considered in Ref. 
are unknown for a dense asymmetric baryon matter it would be meaningless to complicate 
the calculations by the Fermi liquid effects. Therefore we discard residual particle-hole 
interactions and consider the problem in the BCS model. 

The paper is organized as follows. Section II contains some preliminary notes. We 
consider the nonequilibrium gap equations for the case of spin-orbit pairing forces and discuss 
the renormalizations which transform the standard gap equations to a very simple form 
valid near the Fermi surface. In Sec. Ill we consider small deviations of the condensate 
from equilibrium caused by weak external fields and derive the equations for the relevant 
anomalous three-point vertices. In Sec. IV we derive, in the BCS approximation, the 
dispersion equations for eigenmodes of the order parameter and find the whole spectrum of 
the undamped oscillations. Section V represents a discussion of the obtained results and the 
conclusion. Some intermediate and additional calculations are contained in three appendices 
for a more deep understanding. 



Throughout this paper, we use the standard model of weak interactions, the system of 
units h — c—1 and the Boltzmann constant ks — 1. 



II. FORMALISM 

We employ the Matsubara calculation technique for the system in thermal equilibrium and 



and p R = (2k + 1)tvT with k = 0,±1,±2,... is the fermionic Matsubara frequency which 
depends on the temperature T. 

The triplet order parameter, D = D a p (n) , in the neutron superfluid represents a sym- 
metric matrix in spin space (a,f3 =t)i) which can be written as D (n) = Ahag, where 
<r = &2, 03) are Pauli spin matrices, and g = i&2- The angular dependence of the order 
parameter is represented by Cartesian components of the unit vector n = p/p which involves 
the polar angles on the Fermi surface, 



The temperature-dependent gap amplitude A (T) is a real constant (on the Fermi surface), 
and b (n) is a real vector in spin space which we normalize by the condition 



use the standard notation for ordinary propagators of a quasiparticle and a hole, G (p K .p) = 
— *— , G~ (p K .p) = G(—p K ,—p) = — « — , and for anomalous propagators = -* — «- , 
= « » in the momentum representation, where p is the quasiparticle momentum, 



n\ = sin 6 cos ip, ri2 = sin 9 sin ip, n 3 = cos 6. 



(1) 



<P(n)) = l. 



(2) 



Hereafter we use the angle brackets to denote angle averages, 




(3) 



The analytic form of the quasiparticle propagators can be written as 



G (p K , p) = G [p K , p) i, = F (p K , p) bag, 



(4) 



(-Pk, -p) = gv^F {p K , p) , 



(5) 



where we define the scalar Green functions 



The quasiparticle energy is given by 

E p = yje* + AW (n) (7) 

with 

e p ~v F (p- p F ) , (8) 

where Vp <C 1 is the Fermi velocity of the nonrelativistic neutrons. 

At equilibrium, the gap matrix D is connected to the anomalous Green's function F by 
the gap equation [32] 

Da/3 (P) = ~T 22 / ^-T r «7,/35 (P> P') F lS (Pk-P') , (9) 

where T a p tl s (p, p') stands for the block of the interaction diagrams irreducible in the channel 
of two quasiparticles. 

We consider small periodic departure from equilibrium of the form oc exp (zqr — icot). 
By assuming that the interactions between quasiparticles are essentially instantaneous on 
the time scale of the oscillation frequencies the nonequilibrium (time-dependent) gap matri- 
ces Q a f3 (cu, q; p) exp (zqr — iut) must obey the self-consistency conditions (gap equations) 
which depend now on the energy u and space momentum q of the perturbation. Then the 
order parameters D a p (cj, q; p) are to be found as the analytic continuation of the func- 
tions D A (lo v , q; p) connected to the anomalous Green functions J 7 ^' (lu v , q; p K , p) by the 
equations 

-^T ai ,fss (p, p') J^s (^vi pO > ( 10 ) 

re 

s$ (u v , q; p) = —T ^2 j l^/^ (p'> p) K q;p«, p) > (n) 

re 

where a;,, = 2m7]T with 77 = 0, ±1, ±2, ... is the bosonic Matsubara frequency. 

In these equations, the integration goes over infinite momentum space while the quasi- 
particle approach is valid only near the Fermi surface. To get rid of the integration over the 
regions far from the Fermi surface in Eqs. f TTUT) and (1111) . we renormalize the interaction as 



suggested in Refs. [33|], [34|. We define 

r af,l38 (P, P'; T) = r a 7 ,/3<5 (P> P') 

- T E / ^^ ^ ^(p ; p ,, ) (GG-);r; 7 ^(p",p';T), (12) 
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where g = pfM*/tt 2 , the effective mass of a neutron quasiparticle is defined as M* = 
Pf/vf, and the product (GG~)" N = Gn (p K , p") Gn (— p K , — p") is evaluated in the normal 
(nonsuperfluid) state by assuming Gn {p K , p) = (Wk — £ p)~ ■ 

By acting from the left onto the both sides of Eqs. ffTUj) and ffTTl) by the operator 

!- T E/ ^^ U (P,P ^(GG-)'^.., (13) 

and making use of the fact that 



K lt 0S (P. P'; T ) = F « 7 ,/35 (P> P') 

" T E / (P' T ) (^)jv r /n,^ (P"> P) , (14) 



we can recast Eqs. ( TTTTh and (TTTI) to the form 

®S (P) = " T E / S r "7,^ (P. P') - , (15) 

K 

(P) = " T E / (P'> P) - S$GWG>r] p? . (16) 

For brevity we omit the dependence of functions on u v and q. The product GnG n is to be 
evaluated for u v = 0, q = 0. In these equations the integrand decreases very rapidly with 
the distance from the Fermi surface. Since we are interested in the processes near the Fermi 
surface, we can replace all smooth functions of the momentum p with their expressions at 
P = Pf- 

The pairing in the high-density neutron matter is mostly caused by the attractive spin- 
orbit interactions. In the vector notation one can write the renormalized 3 P 2 interaction in 
the form 

2 

QK-iW (p, p) = -V Yl t b ^( n ) • d-sHap l9* ■ Ki( n %s » ( 17 ) 

M=-2 

where V is the constant amplitude taken at the Fermi surface, and bjvf (n) are the vectors 
in spin space which generate standard spin-angle matrices according to 

±=b M (n)*g = (M a/3|1Ms ) ( llM s M ^M) Y 1>Ml (n) . (18) 



M S +M L =M 

These are given by 



b = (-ri!, -n 2 , 2n 3 ) , bi = - y/S/Z (n 3 , in 3 , m + in 2 ) , 

b 2 = (ri! + in 2 , m x - n 2 , 0) , b_ A/ = (-) M b* M . (19) 
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The vectors are mutually orthogonal and are normalized by the condition 



(b* M ,b M ) = 5 



M'M- 



(20) 



It is convenient to divide the integration over the momentum space into integration over 
the solid angle and over the energy according to 



and denote 



d 3 p m . . f dn r°° , 1 
?Tr ... =« / — / de v - Tr 



/oo 
. . -oo 



(21) 



p,K 



Then instead of Eqs. (fl5|) and ([161) we obtain 



£>« (n) = V^b M (n)^ 



M 



X 



frb*) [P {1) -& 1] G N G 



(22) 



^ (n) = V^<xb M (n) 



^-Tr[(b^)^(^ 2 )-S)( 2 )G^ 



(23) 



P,K 



where the dependence of functions on u v and q is again omitted. 

In the absence of external fields the equilibrium order parameters are of the simple form 



Z) (1) (n) = t) (n) = Abag, 

D (2) (n) = D + (-n) = A^o-b. 
From Eqs. ((22) and we obtain 



(24) 
(25) 



b M (n) / b^b ^2(F- AG n G n ) \ . (26) 

M \ p,K I n , 

Using also the fact that F = A (GG~ + b 2 FF) , Eq. (I2"6"j) can be written in the form 

b(n) = V^b M (n)<bb; f A>, (27) 



where the function A is defined as 

/oo 
de p T ( GG ~ ~ G nG n + b 2 FF) . (28) 
-oo „ 

Its explicit form is given in Eq. (|43p . 
With the aid of the expansion 

b(n) = ^<bb; f )b M (n), (29) 

M 

and of Eqs. ([2]) and ([2"01 one can recast Eq. (|2~TI) as 

1 = V (b 2 A) . (30) 

III. VERTEX EQUATIONS 

We now consider the order parameter of the superfluid Fermi-liquid in weak external 
fields {oj, q) exp (iqr — iut). The field interaction with a superfluid should be described 
with the aid of four effective three-point vertices. There are two ordinary effective ver- 
tices, f/j, (uj, q;p) and t~ (p) = fj (— p), corresponding to the creation of a particle and a 
hole by the field (These differ by direction of fermion lines), and two anomalous vertices 
Tji (u>, q; p) corresponding to the creation of two particles or two holes. The latter de- 
scribe the linear departure from equilibrium of the matrix order parameters caused by the 
external field of wave vector q and frequency uj: 

SI)™ (w, q; p) = efW ( Wj q; p ) ^ (co, q) . (31) 

Here the coupling constant e is specified by the external field and a summation over repeated 
Dirac indices \l — 0, 1, 2, 3 is implied. The linearization of Eqs. (122]) and (|23|) by taking into 
account Eq. fl26l) results in equations for the anomalous three-point vertices ' (to, q; n) 

M 

\ PA 

M 



g&b 



M 



5_PV 
eAv 



fPG N G- t 



(32) 



\ eA» 



f^G N G 



(33) 
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where the product GnG^ is to be evaluated for u = 0, q = 0. The corrections to the 
anomalous propagators <5P ( - 1 ' 2 - ) arise owing to the variation of the order parameters and the 
ordinary self-energy of a quasiparticle in the external field and can be depicted symbolically 
by the graphs in Fig. [TJ 




FIG. 1. Graphs corresponding to linear corrections to the anomalous propagators. The external 
field is shown by dashed lines. The anomalous vertices are shown by filled triangles. Filled circles 
are dressed ordinary vertices. 



Analytically one has 



G (P + ) fjpd- (P_) + P (1) (P+) if )p« (P_ 



+ G (P + ) (P_) + P« (P + ) t-G- (P_) , (34) 



and 



^ ~ G~ (P + ) f f G (P_) + P (2) (P + ) T«p( 2 ) (P_) 

+ P( 2 > (P + ) f M G (P_) + G" (P + ) f M P( 2 ) (P_) , (35) 

where we denote P + = {p K + u;^, p + |q) and P_ = p — |q) . 

We shall examine collective excitations of the condensate below the pair-breaking thresh- 
old. The eigenmodes are known to manifest themselves as sharp resonances of the linear 
response of the medium onto external perturbations. Throughout this paper we consider 
the case of 3 P2 neutron pairing, when quasiparticles pair in the most attractive channel with 
spin, orbital and total angular momenta, S — 1, L — 1, J = 2, respectively. This allows one 
to search for the anomalous vertices near the Fermi surface in the form of expansions over 
the eigen functions of the total angular momentum (J, M) with J = 2 and M = 0, ±1, ±2: 

T« {u, q; n) = ^ q) (*b M ) & (36) 

M 



f V ( W| q; n) = 4? (w, q) $ (*b M ) • (37) 



Hereafter we omit the Dirac indices by assuming that the equations are valid for each of 
the vertex components. Inserting expressions (136]) and (137]) into Eqs. f )32|) and (133]) by 
making use of Eq. (1A6[) . one can obtain the set of equations for Bfy ■ The amplitude of 
the pairing interaction V can be eliminated by virtue of the gap equation f )30|) . Substantial 
simplifications can be further done by making use of Eqs. ( 1A3I) . (1A4I) . and ( IA5I) . 
In this way we obtain two uncoupled sets of equations, 

£ 2 < [(ft 2 - r? 2 - 6 2 ) b* M b M , + (b^b) (b M 'b)] X FF ) B m , 

M' 

+ ((b 2 M -b 2 )A)B M = R- M , (38) 



and 



2 < [(^ - r/ 2 ) b* M b M , - (b* M b) (b M ,b)] Iff) B+, 

M> 

+ {{b 2 M -b 2 )A)B + M = R + Ml (39) 



for the linear combinations of the unknown functions 



Bb = l(B$±B$), (40) 

where 

n = 4-^ = ^. ( 41 ) 

6^- = b^ f bjv/ and v is a vector with the magnitude of the Fermi velocity vf and the direction 
of n. The right-hand sides of Eqs. ( 138]) and ( 139]) are given by 

Rm = t\ (Tr [(&h* M ) (fab =F <rbf )] (ft + 77) Z FF ) 

± l - <Tr [(<xb^) (ffb^f-p =F gf-gab)} (ft - r?)Z Fi? ) . (42) 

The function A (n, T) is of the form 

A(n)= [ de ( — tanh — - — tanh — ) , (43) 

v ; \2E 2T 2e 2T J ' v ; 

and the function Xff (w, q; n, T) is defined in Appendix A. For later use we indicate the 

explicit form of this function at q — > 0: 

/DC ^2 

~p T^i — 1 — ^2 tanh Wf, > ( 44 ) 
.oo £/ 4£^ - (a; + «0) ^ 



where 

E = ^e 2 + A 2 b 2 (n) . (45) 

It is necessary to stress that in obtaining Eqs. ( J39l) and (1381) the particular form of the 
gap structure was used, where A (T) is a real scalar and b (n) is a real vector. Therefore 



where 



the present theory does not allow comparison with the theory of 3 He-A by Wolfe 

b (n) = (ni + in 2 ) w with w being a unit vector in spin space. However, the equations can 

_ 

be easily generalized to the case of 3 He-B [30|] where b (n) = n (See Appendix B). 

IV. EIGENMODES OF THE ORDER PARAMETER 

We now look for eigensolutions of Eqs.( l38j) and fl39|) in the BCS approximation, i.e. 
regarding the terms R^ involving f = (f , f) as sources which they are in the limit of 
vanishing Fermi-liquid interactions. We consider the vertices for density, current, spin- 
density and spin-current perturbations interacting with the corresponding vector or axial- 
vector field A M . 

Expansions (1361) and (1371) imply both the unitary and nonunitary excited states. The 
unitary excitations have no spin polarization. In a nonunitary state the Cooper pairs at 



point n have a net average spin 35j, 36]. Evidently the unitary perturbations can be 
excited by the vector external fields while the nonunitary perturbations interact with the 
axial-vector external fields. In the BCS approximation the corresponding ordinary vector 
and axial-vector vertices are of the form f M = (1, v) 1 and (v-<x, <x), respectively. 



A. Unitary excitations 

From Eq. ()42l) we obtain for the vector field. 



R M = (ijL (bb* M l FF ) , 1 <v (qv) (bby l FF )^j , R + M = 0. 



(46) 



In this case from Eq. fl39|) we obtain the trivial solution = 0, and inspection of Eq. 
(|38|) reveals that the eigenmodes which can be excited by the vector external field satisfy 
the dispersion equation 

det |2 < [(n 2 -r, 2 - b 2 ) h* M h M > + (b^b) (b^b)] X FF ) 
+5 M M>({b 2 M -b 2 )A)\=0. (47) 
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We examine the eigenmodes with q = 0. For the 3 P2 pairing with M = one has b = bo- 
Then 

b 2 = 1(1 + Snl) , (48) 

and the functions A (n) = A (n 3 , y) and X (to, n; y) = I (u,n 3 ; y) are axially symmetric at 
q = 0. Here and below we denote y = A (T) /T. In the case q = the azimuth-angle 
integrals can be performed in Eq. fj47|) making use of the orthogonality relations 



dp,* , x , 2 

7T - b M b M' — OmM'Vmi 
Z7T 



(49) 



/dtp 
— (b^b ) (b M /b ) = 8 m ,m> (K/bo) (b M b ) , (50) 

which can be easily verified. As a result after some algebraic manipulations, the dispersion 
equation ( )4"T|) at q = can be obtained in the form 



Y[f M (u,y) = o 



M 



with 



f M (Lu,y) = ((b 2 M -b 2 )A + 2 



to 



(b* M b) (b M b) 



(51) 



(52) 



We thus obtain separate eigenvalue equations for each value of M = 0, ±1, ±2 in the vector 
channel: 



tO 



2A 2 



2 < \b%l - (b* M b) (b M b)] Jo) - < - b 2 ) A) . 



(53) 



At to > 2Ab the function Xo(u,ns,y) has the imaginary part owing to the breaking and 
formation of Cooper pairs. Therefore the undamped collective modes could exist only at 
< to < 2AyJlj2. 

From Eq. ( |53|) we obtain the eigenvalue equation for a perturbation with M = 0, 



lo; 



/ dn 3 (l + 3n 2 3 ) X (u , n 3 ; y) = 0, 
J o 



yielding 



w (j/) = 0. 



(54) 



(55) 



At M = ±1 we obtain the equation 



dn 3 
1 

3./n 



l (1 + 4) - (3nl + Ant + l) 



2A 2 
i 

dn 3 (l - 3^3) A (n 3 ; y) 



X Q (ui,n 3 ;y) 



(56) 
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which has no real solutions in the domain < co < 2Ayl/2. 

For M = ±2 we find u (q = 0) = u 2 (y) which satisfy the dispersion equation 



4A^ 2 



(1 - n 2 3 ) - - (1 - nl) (1 + 7n^) 



X (uj 2 ,n 3 ;y) 



/ dn 3 (l - 3nl) A(n 3 ;y) . 
Jo 



(57) 



B. Nonunitary excitations 

For the axial-vector field from Eq. (T42"|) we find 

1 



R-M — 0, i? 



M 



' I ^ <,(<l v ) ( b A/ xb ) vIff) , ^ <(b; f xb) X FF ) 



(58) 



In this case in Eq. (I3"8"j) the right-hand side vanishes and we obtain the trivial solution 
= 0. Inspection of Eqs. (|39|) reveals that the eigenmodes of the pseudo-vector current 
satisfy the dispersion equation 



det |2 < [(n 2 - V 2 ) b* M b M , - (b* M b) (b M 'b)] Iff) 
+5 MM ,((b 2 M -b 2 )A)\=0. 



(59) 



In the limit q = the manipulations with making use of Eqs. ( H9j) and ( 1501) allow one to 
obtain Eq. fl59|) in the form 

Y[f M (u>,y) = o, (60) 



M 



where 



f M (co,y) = ((b 2 M -b 2 )A + 2 



4A2 



b 2 M - (b^b) (b M b) 



(61) 



The eigenvalues % for M = 0, ±1, ±2 in the axial- vector channel can be found from the 
dispersion equation 

,2 



b M \ 2 l ) = 2 <(b* M b) (b M b) Jo) - <(|M 2 " b 2 ) A) . 



2A 



(62) 



For M = we obtain 



2A 



-a/ 



(i + 3n 2)_(i + 3n 2y 



X (w, y) = 0. 



(63) 



This equation has no real solutions at |u;| < 2A a/1/2. 
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For M = ±1 Eq. (1621) yields (q = 0) = u\ (y) satisfying 



2A 2 



Cj\ (i + nj) - n| (1 - nl) 



X (wi,n 3 ;y) 



dn 3 (l - 3723) A (n 3 ; y) • 



At M = ±2 we find u (q = 0) = U2 (y) satisfying the dispersion equation 



(64) 



dm 



—tf 2 (l-nl)-{l-niy 
1 

dn 3 (l - 3nl) A (n 3 ; y) . 



(a} 2 ,n 3 ;y) 



(65) 



The eigen frequencies of the collective excitations in the limit q — >■ are represented in 
Fig. [2] versus reduced temperature T/T c . In this plot, the oscillation frequencies are shown 
both in units of A = A (T = 0) = 1. 68T C (See Appendix C) and in units of A (T). We 
assume 



13| y (T/Tc) = V2 y/l - T/T c (0.7893 + 1.188T C /T). 



V. DISCUSSION AND CONCLUSION 

Since the equilibrium state of the condensate corresponds to M = the variables Bq ' are 
associated with the broken gauge symmetry. In Fig. [2} the lowest curve uq = represents the 
corresponding Goldstone's mode u> = u (q) at q = 0. The natural appearance of this mode 
is caused by spontaneous breaking of the baryon number owing to the Cooper condensation. 
The collective motion of the condensate in this wave is a periodic variation of the total 
phase without a change of the order parameter structure. The previous analysis [15j has 
shown that, at q > 0, the Goldstone's mode is anisotropic and strongly renormalized by the 
Fermi-liquid interactions. The sound velocity is found strongly dependent on the Landau 
parameters describing the residual particle-hole interactions. 



It was expected 37| that the spontaneous breaking of rotational invariance owing to the 
3 P2 pairing leads to the appearance of three more Goldstone bosons (angulons). We found 
no Goldstone modes associated with the breaking of rotational symmetry. In our model, 
the variables B±{ are associated with a flapping motion of the total angular momentum 
which represents the nonunitary state, where the excited Cooper pairs have a nonzero av- 
erage spin. In other words this can be imagined as a departure of the symmetry axis of the 
bound pair about the symmetry axis of the equilibrated condensate. This is equivalent to 
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FIG. 2. Collective frequencies of the neutron S P 2 condensate at q = in units of Ao = A (T = 0) 
and in units of A (T) versus reduced temperature in the BCS approximation. 

oscillations of the preferred direction of the Cooper pair about the axis of the energy gap in 
the quasiparticle energy. The quasiparticle readjustment can not follow this rapid motion. 
As a result instead of the angulons we found the modes G)\ which correspond to M = ±1 



31] . The temperature variation of uj\ is 



and look like the "normal-flapping" mode in 3 He-A 
naturally explained: the moment of inertia associated with the quasiparticles declines owing 
to a decrease of number of thermal quasiparticles along with a lowering of the temperature. 
At zero temperature the restoring force is expected to vanish along with the number of ther- 
mal quasiparticles. Accordingly, the frequency of these modes passes through a maximum 
and tends to zero for T — > 0. 

The variables B±2 are found as q — > to oscillate at frequencies u = u 2 (y) and u = 
u 2 (y) with u 2 (T = 0) = 1.20A (0) and u 2 (T = 0) = 0.61A (0). To understand the nature of 
these excitations we introduce three spontaneous orbital axes of the gap j, k, 1 with j _L k and 
1 = j x k being the symmetry axis. Then the equilibrium gap matrix of the 3 P 2 condensate 
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with M = at the point n on the Fermi surface can be written as 

(j — ik) n 21n 
21n - ( j + ik.) n 

Variation of this matrix caused by the modes with M = ±2 is of the form (we omit the 
superscript in B\) 




8D\ 2 \ = 



-B 2 (j + ik) n 

£L 2 (j - ik) n 

This form indicates that the oscillations occur in the plane (j,k). The admixture of the 
5D 2 to D can be imagined as a clapping motion of j and k about their equilibrium relative 
angle of n/2. 

According to Eqs. (HB"]) . in the vector channel of excitation one has B^ 2 = B 2 , and the 
excited state D — D + 6D\ 2 \ is unitary. This means the u 2 excitations correspond to the 
states with no net spin polarization. In the axial-vector channel, from Eqs. (I58j) one obtains 
£>_ 2 = —B 2 , and the excited state D = D + 5D\ 2 \ is nonunitary, i.e. for u 2 excitations the 
average spin expectation value is nonzero. (The same takes place for the flapping mode cDi.) 
This allows one to term these oscillations by spin waves. 

The flapping mode u\ needs special discussion. The oscillation frequency for this mode is 
small, uf <C 4A 2 . This allows one to neglect co 2 in the integrand of Eq. (jUJ) thus obtaining 

^ /- \ t / \ A 2 f°° de E 

T (wi, n) ~ J (n) = — / — tanh — (66) 



with E = ^e 2 + A 2 b 2 (n 3 ). 

After this simplification we find the analytic solution of the form 



2 / dn 3 
0? = A 2 -l? 



3n 2 (l-n 2 )J -(l-3n|)A 



(67) 

3 Jo dn 3 ( 1 + n i) I o 

In previous works [17J the right-hand side of Eq. fl67|) was evaluated in the average-angle 
approximation assuming that the anisotropic gap in the quasiparticle energy is replaced by 
its average-angle magnitude, 6 2 A 2 (T) — > A 2 (T). Such approach results in = A (T) /y5 
with a simple temperature dependence of the excitation energy only through the gap ampli- 
tude. Owing to the gap anisotropy the energy of the spin density oscillations depends more 
dramatically on the temperature, as shown in Fig. |2j As can be seen the gap anisotropy 
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leads to a strong decreasing of the level energy along with a lowering of the temperature. 
This is to suppress substantially the neutrino energy losses caused by the spin wave decays 
because the rate of neutrino losses is strongly dependent on the wave energy. Since we have 
found also the new collective modes, which are kinematically able to decay into neutrino 
pairs, the problem of neutrino emissivity of the 3 P2 condensate at lowest temperatures is to 
be revisited. 

It is necessary to note finally that the Fermi-liquid effects can somewhat modify the 
collective frequencies as compared to that found in the BCS approximation. 

Appendix A: Functions used 

We denote as Lx,x (w, q; p) the analytical continuation of the Matsubara sums: 

W ("„ P+f ; P-f ) = T J2 x (P+) X' (P_) , (Al) 

K 

where X, X' G G, F, G~ and operate with integrals over the quasiparticle energy: 

X XX ' (w, n, q; T) = J de p L X x' (u, p+^, P~|) ■ ( A2 ) 

These are functions of u, q and the direction of a quasiparticle momentum n, except the 
function X GnG - which is to be evaluated for u = 0, q = in all cases. One can easily check 
that 

%G-G — %GG- j %GF — —%FG > ^G~F = ~^FG- ■ (A3) 

We use also the following relations valid for arbitrary u, q, T 

1g-f = ^^Zff , (A4) 
u + qv 

FG = — 2A — ' ^ ^ 

and 

Zgg- ~ I GnG - + b 2 l FF = A+ U " ~^ V) 1 FF , (A6) 
The function A (n) is defined by Eq. (1251) and can be written as 



Zgg- ~ x g n g n + b2 lFF 



A(n) 

The explicit form of this function is given in Eq. (TJ3 



w =0,q=0 



(A7) 
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Appendix B: Eigenmodes in 3 He-B 



In the liquid 3 He-B, the pairing occurs owing to the exchange central interaction which is 
independent of the total angular momentum of the bound pair, J. In this case the interaction 
in the channel S — 1, L — 1 can be written in the same form as given in Eq. (|17l) but the 
summation over the total angular momentum J = 0, 1, 2 is to be added with — J < M < J. 
Then instead of Eqs. (|47|) and (|59|) one obtains, respectively, 



and 



det | < (bj/ M 'b} M - Sjj^MM'h 2 ) A) 



+ 



to 2 - (qv) 2 



2A 2 



2b 2 (Joj, M ,b* JM ) + 2 (b* JM b) (b^^b) 



X 



(Bl) 



det | ( (bj/ M 'b} M - Sjj^MM'b 2 ) A) 



+ 



u 2 - (qv) 2 
2A2 



( b J'A/'b} M ) - 2 (b} M b) (bj/ M /b) 



where the vectors hjM for J = 2 are given by Eq. (fl9l) and for J = 0, 1 one has 



(B2) 



b ,o = (-ni, -n 2 , -n 3 ) , 
^3 



b i,o = \l^( in 2, -inx,0) 



bi,i = b* _! = ^/- (-n 3 , -m 3 , m + in 2 ) , 



(B3) 



In 3 He-B, the equilibrium state of the condensate corresponds to the bound pairs with 
J = 0. The order parameter in such a system is h = bo,o, and the energy gap Ab (n) = A 
is isotropic (b 2 = 1). The latter means that, in Eqs. (IBlj) . (1B2j) . the functions A and X are 
isotropic and can be moved beyond the integrals. Using also the orthogonality condition, 
(b JM b}, M ,) = Sjj'Smw, we find 



J'M'Ujm ~ °JJ'°MM 



S.jj'SMM'b ) = 



(B4) 



thus obtaining from Eq. ( 1B1I) the dispersion equation 



det 



2A 2 



2 (bj, M ,b} M ) + 2 (b} M b) (bj/^/b) 



0. 



(B5) 
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We examine solutions to this equation for q = 0. A straightforward calculation of the angle 
integrals results in 

(B6) 



(B7) 



((b} M b) (bj/M'b)) = Sjj'Smw [Sjo^mo + ^.n 
Implying J = 0, 1, 2 and \M\ < J we obtain the dispersion equation 

Q 2 (fi 2 -l) 3 ^ 2 -0 5 = O , 



where we denote 



n 



oo 



2A(T)- (B8) 
From this equation we find one mode with frequency oo = 0, five degenerate modes with 



frequency oo = \J 12/5A, and three degenerate modes with frequency oo = 2A. In the same 
way starting from Eq. (IB 2 1) we obtain the dispersion equation for eigenmodes of spin density 
and spin-current density which can be excited by the axial- vector external field 



(^ 2 -i) (^) 3 (^ 2 -0 5 = O- 



(B9) 



From this equation we find one mode with frequency oo = 2 A, five degenerate modes with 
frequency oo = a/8/5A, and three degenerate modes with frequency oo = 0. This completes 
the enumeration of the 18 B-phase collective modes for q = 0. Note these results are identical 
to ones obtained earlier in Refs. 30], 36 ]. 



Appendix C: Critical temperature 

Combining Eq. ( 130]) for T > with the same equation for T = one can obtain the 
equality 

(, °° rlr i N 

Wfc / T?. IT . »\— 1 



de 



{e 2 + A 2 b 2 y 1/2 - {e 2 + A 2 b 2 ) 



2Z2\-l/2 



i(^ /r +iy 



(Cl) 



where E = V £ 2 + A 2 b 2 . Performing the integral in the left-hand side we find 



ln^ = 2( !r 



dx 



A \ Jo a/x 2 + y 2 b 2 
The internal integral in the right-hand side is calculated in Ref. 



point T — > T c one obtains 



A o /r 2 A„ *T , 7C(3)A 2 



ln X = \ 6 i ln ^Al + Svr 2 T 2 



(C2) 

Close to the transition 
(C3) 
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where In 7 = C = 0.577 is Euler's constant. According to this equation the gap vanishes at 
the temperature 

T c = A ^exp(6 2 ln6). (C4) 

71 

We evaluate this relation for b 2 = (1 + 3nl) /2 to find 



A 7T 



e27^ 3 -2 = 0.595. (C5) 



[1] P. Haensel, Nucl. Phys. A 301, 53 (1978). 

[2] L. B. Leinson, V. N. Oraevsky, and V. B. Semikoz, Physics Lett. B, 209, 80 (1988). 
[3] C. J. Horowitz and K. Weinberger, Nucl. Phys. A 531, 665 (1991). 
[4] C. J. Horowitz and K. Weinberger, Phys. Lett. B 266, 236 (1991). 
[5] E. Braaten and D. Segel, Phys.Rev. D 48, 1478 (1993). 
[6] F. Matera and V. Yu. Denisov, Phys. Rev. C 49, 2816 (1994). 
[7] S. Reddy, M. Prakash, J. M. Lattimer, and J. Pons, Phys. Rev. C 59, 2888 (1999). 
[8] L. B. Leinson, Phys. Lett. B 469, 166 (1999). 
[9] L. B. Leinson, Phys. Lett. B 473, 318 (2000). 
[10] L. B. Leinson, Nucl. Phys. A 687, 489 (2001). 

[11] V. Greco, M. Colonna, M. Di Toro, and F. Matera, Phys. Rev. C 67, 015203 (2003). 
[12] C. Providencia, L. Brito, A. M. S. Santos, D. P. Menezes, and S. S. Avancini, Phys. Rev. C. 
74, 045802 (2006). 

[13] D. G. Yakovlev, A. D. Kaminker, and K. P. Levenfish, Astron. Astrophys. 343, 650 (1999). 

[14] L. B. Leinson, Phys. Rev. C 81, 025501 (2010). 

[15] L. B. Leinson, Phys. Rev. C 83, 055803 (2011). 

[16] L. B. Leinson, Phys. Lett. B 689, 60 (2010). 

[17] L. B. Leinson, Phys. Rev. C 82, 065503 (2010). 

[18] L. B. Leinson, Phys. Lett. B 702, 422 (2011). 

[19] L. B. Leinson, Phys. Rev. C 84, 045501 (2011). 

[20] M. Hoffberg, A. E. Glassgold, R. W. Richardson, M. Ruderman, Phys. Rev. Lett. 24, 775 
(1970). 

19 



M.V. Zverev, J. W. Clark, and V. A. Khodel, Nucl. Phys. A 720, 20 (2003). 

V. A. Khodel, J. W. Clark, and M.V. Zverev, Phys. Rev. Lett. 87, 031103 (2001). 



V. A. Khodel, J. W. Clark, and M.V. Zverev, e:Print: |arXiv:nucl- th/0203046 , 
P. Wolfe, Phys. Rev. Lett. 30, 1169 (1973). 
P. Wolfe, Phys. Rev. Lett. 31, 1437 (1973). 
K. Maki and H. Ebisawa, J. Low Temp. Phys. 15, 213 (1974). 
R. Combescot, Phys. Rev. A 10, 1700 (1974). 
R. Combescot, Phys. Rev. Lett. 33, 946 (1974). 
H. Ebisawa and K. Maki, Prog. Theor. Phys. 51, 337 (1974). 
P. Wolfe, Physica B 90, 96 (1977). 
P. Wolfe, Phys. Rev. Lett. 37, 1279 (1976). 
L. Amundsen and E. 0stgaard, Nucl. Phys. A 442, 163 (1985). 

A. I. Larkin and A. B. Migdal, Zh. Eksp. Teor. Fiz. 44, 1703 (1963) [Sov. Phys. JETP 17, 
1146 (1963)]. 

A. J. Leggett, Phys. Rev. 147, 119 (1966). 
A. J. Leggett, Rev. Mod. Phys, 47, 331 (1975) 

J. B. Ketterson and S. N. Song, Superconductivity, (University press, Cambridge, 1999), p. 
414. 

P. F. Bedaque, G. Rupak, and M. J. Savage, Phys. Rev. C 68, 065802 (2003). 

E. M. Lifshitz, L. P. Pitaevskii, Statistical Physics, Part 2 (Pergamon Press, Oxford, 1980). 



20 



